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Abstract. We review the most recent evidence for the amazing properties of the 
density distribution of the dark matter around spiral galaxies. Their rotation curves, 
coadded according to the galaxy luminosity, conform to an Universal profile which can 
be represented as the sum of an exponential thin disk term plus a spherical halo term 
with a flat density core. From dwarfs to giants, these halos feature a constant density 
region of size ro and core density po related by po — 4.5 x 10~'^(ro/kpc)~^'''^MQpc~''. 
At the highest masses po decreases exponentially with ro, revealing a lack of objects 
with disk masses > lO^^M© and central densities > 1.5 x 10~^ (ro/kpc)~^M0pc~^ 
implying a maximum mass of « 2 x for a dark halo hosting a stellar disk. 

The fine structure of dark matter halos is obtained from the kinematics of a number of 
suitable low-luminosity disk galaxies. The halo circular velocity increases linearly with 
radius out to the edge of the stellar disk, implying a constant dark halo density over the 
entire disk region. The properties of halos around normal spirals provide substantial 
evidence of a discrepancy between the mass distributions predicted in the Cold Dark 
Matter scenario and those actually detected around galaxies. 



1 Introduction 

Rotation curves (RC's) of disk galaxies are the best probe for dark matter (DM) 
on galactic scale. Notvifithstanding the impressive amount of knowledge gathered 
in the past 20 years, only very recently we start to shed light to crucial aspects 
of the mass distribution of dark halos, including their radial density profile, 
and its claimed universality. On a cosmological side, high-resolution N-body 
simulations have shown that cold dark matter (CDM) halos achieve a specific 
equilibrium density profile [13 hereafter NFW, 5, 8, 12, 9] characterized by one 
free parameter, e.g. the halo mass. In the inner region the DM halos density 
profiles show some scatter around an average profile which is characterized by a 
power-law cusp p ~ r~^, with 7 = 1 — 1.5 [13, 12, 2]. In detail, the DM density 
profile is: 

PNFw{r) = — — ^ (1) 

where rj and ps are respectively the characteristic inner radius and density. Let 
us define rvir as the radius within which the mean density is Z\vir times the 
mean universal density pm at the halo formation rcdshift, and the associated 
virial mass Afvir and velocity Vvw = GMvir/rvW- Hereafter we assume the ylCDM 
scenario, with fim — 0.3, ]7/i — 0.7 and h — 0.75, so that Zivir — 340 at z ~ 0. By 
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Fig. 1. Synthetic rotation curves (filled circles with error bars) and URC (solid line) 
with the separate dark/luminous contributions (dotted line: disks; dashed line: halos). 



assuming the concentration parameter as c = r^ir/rs the halo circular velocity 
^Fw('') takes the form [2]: 

^NFW — ^vir '^^J^ ^ ^ (2) 

where x = r/vg and A{x) = ln(l + x) — x/{l + x). As the relation between 
Vyir and Tvir is fully specified by the background cosmology, the independent 
parameters characterizing the model reduce from three to two (c and rg). Let 
us stress that a high density fim = 1 model, with a concentration parameter 
c > 12, is definitely unable to account for the observed galaxy kinematics [11]. 

So far, due to the limited number of suitable RC's and to the serious un- 
certainties in deriving the actual amount of luminous matter inside the inner 
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regions of spirals, it has been difficult to investigate the internal structure of 

dark hales. These difficulties have been overcome by means of: 

i) a specific investigation of the Universal Rotation Curve [16], built by coadding 
1000 RC's, in which we adopt a general halo mass distribution: 

V^URci^) = Kit (1 - /3) (1 + «') (^2 + ^2) (3) 

with X = r/ropt, a the halo core radius in units of rgpt and /? = (Vafj//jc'(''opt)/Kpt)^ 
It is important to remark that, out to rgpt, this mass model is neutral with re- 
spect to the halo profile. Indeed, by varying /3 and a, we can efficiently reproduce 
the maximum disk, the solid -body, the no-halo, the all halo, the CDM and the 
core-less-halo models. For instance, CDM halos with concentration parameter 
c = 5 and = ropt are well fit by (3) with a ~ 0.33 

ii) a number of suitably selected individual RC's [1], whose mass decomposition 
has been made adopting the cored Burkert-Borriello-Salucci (BBS) halo profile 
(see below). 

2 Dark Matter Properties from the Universal Rotation 
Curve 

The observational framework is the following: a) the mass in spirals is distributed 
according to the Inner Baryon Dominance (IBD) regime [16] : there is a character- 
istic transition radius rjBD — 2r£i(V^pt/220 km/s)^'^ (r^ is the disk scale-length 
and Vopt = V{ropt)) according which, for r < riBD, the luminous matter to- 
tally accounts for the gravitating mass, whereas, for r > rjBD, the dark matter 
shows dynamically up and rapidly becomes the dominant component [20, 18, 1]. 
Then, although dark halo might extend down to the galaxy centers, it is only for 
r > Tied that they give a non-negligible contribution to the circular velocity, h) 
DM is distributed in a different way with respect to any of the various baryonic 
components [16, 6], and c) HI contribution to the circular velocity at r < Vopt, 
is negligible [e.g. 17]. 

2.1 Mass modeling 

Persic, Salucci and Stel [16] have derived from ~ 15000 velocity measurements 
of ~ 1000 RC's the synthetic rotation velocities of spirals Vsyn{:;r—^] j^), sorted 
by luminosity (Fig. 1, with Lj the /-band luminosity (L//i* = l0-(*^^+2i-9)/5). 
Remarkably, individual RC's have a negligible variance with respect to their 
corresponding synthetic curves: spirals sweep a very narrow locus in the RC- 
profile/amplitude/luminosity space. In addition, kinematical properties of spirals 
do significantly change with galaxy luminosity [e.g. 16], then it is natural to relate 
their mass distribution with this quantity. The whole set of synthetic RC's define 
the Universal Rotation Curve (URC), composed by the sum of two terms: a) an 
exponential thin disk with circular velocity (see [16]): 



yluRci^) = 1-28 P {I^K^ - hK^)\i.6. 



(4) 
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Fig. 2. a vs. (3 and /3 Vopt- 



and a spherical halo, whose velocity contribution is given by (3). At high lumi- 
nosities, the contribution from a bulge component has also been considered. 

The data (i.e. the synthetic curves Vgyn) select the actual model out of this 
family, by setting V^j^^ix) = V^^^j^^ix, P,a) + Vl^iic{x,P) with a and f3 as 
free parameters. An extremely good fit occurs for a ~ 1.5(L//L*) [16] or, cquiv- 
alently, for a = a{P) and /? = P{log Vopt ) as plotted in Fig. 2. With these values 
the URC reproduces the data Vsyn{r) up to their rms (i.e. within 2%). Moreover, 
at fixed luminosity the a fitting uncertainties in a and f3 arc lesser than 20%. 
The emerging picture is: i) smaller objects have more fractional amount of dark 
matter (inside ropt: M*/Mvir ^ 0.2 (M*/2 x 10^'^ Mq)^-"^^ [20]), ii) dark mass 
increses with radius much more that linearly. 



2.2 Halo Density Profiles 

The above evidence calls for a quite specific DM density profile; we adopt the 
BBS halo mass distribution [3, 4, 1]: 

3 

(r + ro)(H +r5) 

where po and ro are free parameters which represent the central DM density 

and the core radius. Of course, for ro <S r,;, we recover a cuspy profile. Within 
spherical symmetry, the mass distribution is given by: 

MBBs(r-) = 4 Mo {ln(l + r/ro) - arctan(r/ro) + 0.5 ln[l + {r/rof]} (6) 

Mo ~ 1.6 po rl (7) 

with Mq the dark mass within the core. The halo contribution to the circular 
velocity is then: V^Qg{r) = G Mbbs(^)/^- Although the dark matter "core" 
parameters Tq, po and Mq are in principle independent, the observations reveal 
a quite strong correlation among them [e.g. 19]. Then, dark halos may be an 
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Fig. 3. up) Central halo density po (in g/cm'^) vs. disk mass (in solar units) for normal 
spirals [filled circles); bottom) central density vs. core radii (in kpc) for normal spirals. 
The straight lines are from [3], whereas the curved lines are the best fits used in §4. 



1-parameter family, completely specified by e.g. their core mass Afg- When we 
test the disk+BBS velocities with po and left as free parameters, we find that, 
at any luminosity and out to ~ 6 r^;, the model is indistinguishable from data 
(i.e. Vsyn{f))- More specifically, we reproduce the synthetic rotation curves at 
the level of their rms. The values of and po derived in this way agree with 
the extrapolation at high masses of the scaling law p oc r^'^^^ [3] established 
for objects with much smaller core radii and stellar masses (see Fig. 3). Let 
us notice that the core radii are pretty large (rp r^): ever-rising halo RC's 
cannot be excluded by the data. Moreover, spirals lie on the extrapolation of 

— 1/3 

the disk-mass vs. central halo density relationship po oc found for dwarf 

galaxies [3] , to indicate that the densest halos harbor the least massive disks (see 
Fig. 3). 

The curvature in po vs. tq at the highest masses/lowest densities can be linked 
to the existence of an upper mass limit in M^ir which is evident by the sudden 
decline of the baryonic mass function of disk galaxies at M™"^ — 2x 10^^ Mq [20] . 
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In fact, such a limit implies a maximum halo mass of ~ Qq/Qi, M^"^. 

Then, for (6) and (7), Afvh- = r] Mo, with 77 ~ 12 for (r2o, ,2) = (0.3, 0.03, 3), 
and the limiting halo mass implies a lack of objects with po > 4 x 10~^^ g/cvc? 
and ro > 30 kpc, as is evident in Fig. 3. On the other side, the observed deficit 
of objects with Md ^ M™"^ and po > 4 x 10~^^ g/cm^, suggests that, at 
this mass scale, the total-to-baryonic density ratio nears the cosmological value 
Qo/Qb ~ 10. 

2.3 Testing CDM with the URC 

The BBS density profile reproduces in synthetic RC's the DM halo contributions, 
at least out to two optical radii. This is in contradiction with CDM halo prop- 
erties according to which the velocity dispersion a of the dark matter particles 
decreases towards the center to reach ^ for r ^ 0. Dark halos therefore, are 
not kinematically cold structures, but "warm" regions of sizes ro cx ^'^ which, 
by the way, turn up quite large: ro ~ 4 — 7 r^. Then, the boundary of the core 
region is well beyond the region of the stellar disk and there is not evidence that 
dark halos converge to a p ~ r~^ (or a steeper) regime, as dictated by CDM 
predictions. 

3 Dark Matter Properties from Individual Rotation 
Curves 

Although deriving halo densities from individual RC's is certainly complicated, 
the belief according to which one always gets ambiguous halo mass modeling 
[e.g. 22] is incorrect. In fact, this is true only for rotation curves of low spatial 
resolution, i.e. with less than ~ 3 measures per exponential disk length-scale 
occurring in most HI RC's. In this case, since the parameters of the galaxy 
structure are very sensitive to the shape of the rotation curve in the region 
< r < rd, there are no sufficient data to constrain models. 

In the case of high-quality optical RC's tens of independent measurements in 
the critical region make possible to infer the halo mass distribution. Moreover, 
since the dark component can be better traced when the disk contributes to 
the dynamics in a modest way, a convenient strategy leads to investigate DM- 
dominated objects, like dwarf and low surface brightness (LSB) galaxies. It is 
well known that for the latter [e.g. 7, 11, 3, 4, 9, 10, 21] the results are far from 
being definitive in that they are 1 ) affected by a quite low spatial resolution and 
2) uncertain, due to the limited amount of available kinematical data [e.g. 23]. 

Since most of the properties of cosmological halos arc claimed universal, 
an useful strategy is to investigate a number of high-quality optical rotation 
curves of low luminosity late- type spirals, with /-band absolute magnitudes 
-21.4 < Ml < -20.0 and 100 < Kpt < 170 km s~^ Objects in this luminos- 
ity/velocity range are DM dominated [e.g. 20] but their RC's, measured at an 
angular resolution of 2", have an excellent spatial resolution of ^ 100(Z?/10 Mpc) 
pc and ndata ~ ^opt/w independent measurements. For nearby galaxies: w «rd 
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Fig. 4. Halo parameters (o is in units of ropt) for the individual RC's. Notice that 
for this sample we can derive in model-independent way o ^ 1, in disagreement with 

CDM predictions. 

and Uiiata > 25. Moreover, we select RC's of bulge-less systems, so that the stel- 
lar disk is the only baryonic component for r ^ r^. 

In detail, we extract the best 9 rotation curves, from the 'excellent' subsample 
of 80 rotation curves of [15], which are all suitable for an accurate mass modeling. 
In fact, these RC's trace properly the gravitational potential in that: 1) data 
extend at least out to the optical radius, 2) they are smooth and symmetric, 3) 
they have small rms, 4) they have high spatial resolution and a homogeneous 
radial data coverage, i.e. about 30 — 100 data points homogeneously distributed 
with radius and between the two arms. The 9 extracted galaxies arc of low 
luminosity (5 x lO^L© < Lj < 2 x 10^°Lq; 100 < Vopt < 170 km s"^) and 
their 7-band surface luminosity profiles are (almost) perfect radial exponential. 
These two last criteria, not indispensable to perform the mass decomposition, 
help inferring the dark halo density distribution. Each RC has 7—15 velocity 
points inside ropt, each one being the average of 2 — 6 independent data. The 
RC's spatial resolution is better than 1/20 ropi, the velocity rms is about 3% 
and the RC's logarithmic derivative is generally known within about 0.05. 

3.1 Halo Density Profiles 

We model the mass distribution as the sum of two components: a stellar disk 
and a spherical dark halo. By assuming centrifugal equilibrium under the action 
of the gravitational potential, the observed circular velocity can be split into the 
two components: V^{r) = V^{r) + V^{r). By selection, the objects are bulge-less 
and stars are distributed like an exponential thin disk. Light traces the mass via 
an assumed radially constant mass-to-light ratio. 
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Fig. 5. BBS fits {thick solid line) to tlie RC's {points with errorbars). Tliin solid lines 
represent the disk and halo contributions. Notice the steep halo velocity profiles. The 
maximum and minimum disk solutions {dashed lines) provide the theoretical uncer- 
tainties. 

We neglect the gas contribution Vgas{r) since in normal spirals it is usually 
modest within the optical region [17, Fig. 4.13]: flgas = (^gis/^^)i-<,pt ~ 0.1. Fur- 
thermore, high resolution HI observations show that in low luminosity spirals: 
Vgasir) ~ for r < rd and Vgas{r) ^ (20±5)(r-rd)/2rd for ra < r < 3u. Thus, 
in the optical region: i) V^ag{r) « V'^{r) and ii) d[V'^{r)-V^a^{r)]/dr ^ 0. This 
last condition implies that by including Vgas the halo velocity profiles would re- 
sult steeper and then the core radius in the halo density even larger. Incidentally, 
this is not the case for dwarfs and LSBs: most of their kinematics is affected by 
the HI disk gravitational pull in such a way that neglecting it could bias the de- 
termination of the DM density. The circular velocity profile of the disk is given 
by (4) and the DM halo will have the form given by (3). Since we normalize 
(at ropt) the velocity model {V^ + V^)^^^ to the observed rotation speed Vopt, /3 
enters explicitly in the halo velocity model and this reduces the free parameters 
of the mass model to two. 

For each galaxy, wc determine the values of the parameters /3 and a by means 
of a x^-minimization fit to the observed rotation curves: 



VLdeiir; a) = Vi{r; (3) + V^{r; /3, a) 



(8) 
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Fig. 6. NFW best-fits solid lines of the rotation curves (filled circles) compared with 
the BBS fits (dashed lines). The values are also indicated. 



A central role in discriminating among the different mass decompositions is 
played by the derivative of the velocity field dV/dr. It has been shown [e.g. 
14] that by taking into account the logarithmic gradient of the circiilar velocity 
field defined as: V(r) = ^diogr^^ ""^^ retrieve the crucial information stored 
in the shape of the rotation curve. Then, we set the x^'s as the sum of those 
evaluated on velocities and on logarithmic gradients: Xy = X^"=i ^'"^'""^'j^'^'"*'^'"^ 

and Xv = J27=i with Vmodeiiri, f3,a) given from the above 

equations. The parameters of the mass models are finally obtained by minimizing 
the quantity xLt = Xv + Xv- 

The bc;st fit models parameters are shown in Fig. 4. The disk-contribution 
(3 and the halo core radius a span a range from 0.1 to 0.5 and from 0.8 to 2.5, 
respectively. They are pretty well constrained in a small and continuous region of 
the (a, /?) space. The derived mass models are shown in Fig. 5, alongside with the 
separate disk and halo contributions. We also get a "lowest" and a "highest" halo 
velocity curve (dashed lines in figure) by subtracting from V{r) the maximum 
and the minimum disk c;ontributions Vd{r) obtained by substituting in (4) the 
parameter /3 with /3best + and /Sbest — <5/3> respectively. It is obvious that the 
halo curves are steadily increasing, almost linearly, out to the last data point. In 
each object the uniqueness of the resulting halo velocity model can be realized 
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by the fact that the maximum-disk and minimum-disk models almost coincide. 

Remarkably, wc find that the size of the halo density core is always greater than 
the disk characteristic scale-length and it can extend beyond the disk edge 
(and the region investigated). 

3.2 Testing CDM 

In Fig. 5 we show how the halo velocity profiles of the nine galaxies rise almost 
linearly with radius, at least out to the disk edge: 

Vh{r) oc r 0.05 Vopt ^ Vopt (9) 

The halo density profile has a well defined (core) radius within which the density 

is approximately constant. This is inconsistent with the singular halo density 
distribution emerging in the Cold Dark Matter (CDM) halo formation scenario. 
More precisely, since the CDM halos are, at small radii, likely more cuspy than 
the NFW profile: pcDM r~^-^ [e.g. 14], the steepest CDM halo velocity profile 
Vh{r) (X r^/'* results too shallow with respect to observations. Although the mass 
models of (3) converge to a distribution with an inner core rather than with a 
central spike, it is worth, given the importance of such result, checking in a direct 
way the (in) compatibility of the CDM models with galaxy kinematics. 

So, we assume the NFW functional form for the halo density given by (1), 
leaving c and as free independent parameters, although N-body simulations 
and semi-analytic investigations indicate that they correlate. This choice to in- 
crease the chance of a good fit. We also imposed to the object under study 
a conservative halo mass upper limit of 2 x IO^^Mq. The fits to the data are 
shown in Fig. 6, together with the BBS fits: for seven out of nine objects the 
NFW models are unacceptably worse than the BBS solutions. Moreover, in all 
objects, the CDM virial mass is too high: Myir ^ 2 x IO^^Mq and the resulting 
disk mass-to-light ratio too low ( ^ 10^^ in the /-band). 

4 The Intriguing Evidence from Dark Matter Halos 

The dark halos around spirals emerge as an one-parameter family; the order 
parameter (either the central density or the core radius) correlates with the lu- 
minous mass. However, we do not know how it is related to the global structural 
properties of the dark halo, like the virial radius or the virial mass, unless we 
extrapolate out over the BBS profile. That is because the halo RC, out to the 
outermost data, is completely determined by physical parameters, the central 
core density and the core radius, which have little counterpart in the gravita- 
tional instability /hierarchical clustering picture. 

Caveat the above extrapolation, the location of spiral galaxies in the virtual 
space of virial mass/halo "central" density /stellar mass, that, on theorethical 
basis, should be roughly 3-D random determined by several different and non- 
linear physical processes, is remarkably found to degenerate and to lie on a 
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Fig. 7. The dark-to-luminous mass ratio as function of the normalized radius and the 
total disk mass. 

curve. Indeed, in Fig. 7 we show the dark-to-luminous mass ratio as function of 
the normalized radius and the total disk mass. The surface has been obtained 
by adopting the correlations between the halo and the disk parameters we found 
in our previous works (see Fig. 3 [19]): 

logro = 9.10 + 0.28 logpo - 3.49 x 10^° (10) 

logpo = -23.0 - 0.077 \ogMd - 9.98 x 10"^ M^ '*^ (11) 

and: 

logrd = 4.96 - 1.17 logM^ + 0.070 (logMd)^ (12) 

from data in [16]. The dark-to-luminous mass ratio at fixed ratio increases as 
the total disk mass decreases; for example at r = 10 it raises from 20% for 
massive disks {Md = lO^^j^f^) to 220% for smaller disks (M^ = lO^M©). 

Two conclusive statements can be drawn: dark matter halos have an inner 
constant-density region, whose size exceeds the stellar disk length-scale. Second, 
there is no evidence that dark halos converge, at large radii, to a p ^ r^^ (or 
steeper) profile. 

The existence of a region of "constant" density po — ^ and size ro 

is hardly explained within current theories of galaxy formation. Moreover, the 
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evidence of a smooth halo profile is growing more an more in recent literature 
[e.g. 26, 27, 28, 29] and a number of different solutions have been proposed to 
solve this problem [e.g. 30, 31, 32, 33]. Let us stress, however, that we should 
incorporate all the intriguing halo properties described in this review. 
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